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ABSTRACT 

We present new results on the Hubble diagram of distant type la supernovae (SNela) 
segregated according to the type of host galaxy. This makes it possible to check earlier 
evidence for a cosmological constant by explicitly comparing SNe residing in galaxies 
likely to contain negligible dust with the larger sample. The cosmological parameters 
derived from these SNe la hosted by presumed dust-free early-type galaxies supports 
earlier claims for a cosmological constant, which we demonstrate at ~ 5cr significance, 
and the internal extinction implied is small even for late-type systems {Ab < 0.2). 
Thus, our data demonstrate that host galaxy extinction is unlikely to systematically 
dim distant SNe la in a manner that would produce a spurious cosmological constant. 
We classify the host galaxies of 39 distant SNe discovered by the Supernova Cosmology 
Project (SCP) using the combination of Hubble Space Telescope STIS imaging, Keck- 
II ESI spectroscopy and ground-based broad-band photometry. The distant data are 
analysed in comparison with a low-redshift sample of 25 SNela re-calibrated according 
to the precepts of the SCP. The scatter observed in the SNela Hubble diagrams 
correlates closely with host galaxy morphology. We find, as expected, this scatter is 
smallest for SNe la occurring in early-type hosts and largest for those occurring in late- 
type galaxies. Moreover, SNe residing in early-type hosts appear only ~ 0.14±0.09mag 
brighter in their light-curve-width-corrected luminosity than those in late-type hosts, 
implying only a modest amount of dust extinction even in the late-type systems. As in 
previous studies, these results are broadly independent of whether corrections based 
upon SN light-curve shapes are performed. We also use our high-redshift dataset to 
search for morphological dependencies in the SNe light-curves, as are sometimes seen 
in low-redshift samples. No signihcant trends are found, possibly because the range of 
light-curve widths is too limited. 

Key -words: cosmological parameters — distance scale — cosmology: observations 
— supernovae: general 


1 INTRODUCTION 
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of the cosmoloffical world-model. Althoueh their 

further strengthens these conclusions (see e.g. Bahcall et al. 


physical properties are still not fully understood, by ex¬ 
ploiting an empirical relationship between th e peak intrin¬ 


sic luminosity and the light-curve decay time ( Phillips 


1993 


Riess et al. 

199E|, 1996; 

Hamuy et al. 

19961; 

Perlmutter et al. 

1997), the dispersion ir 

their photometric properties can be 


reduced to ~ 0.17 mag, making them valua ble ‘calibrated 


candle s’ over a wide range in redshift (e.g. Hamuy et al 


1996c). 


Although numerous alternativ e probes of the cosmolog¬ 
ical p arameters are now availa ble ( Jaffe et al, ^001 JPeacock 


et al. 2001; ie Bernardis et al, 2002; [Efstathiou et al. 2002)^ 


in principle the SNIa Hubble diagram provides the only 
direct measure of the cosmic expansion history free from 
any model a ssumptions. The systematic detection (Perl- 


mutter 


1997) and photometric calibration of high redshift 


SNela has led to two independent SNIa Hubble diagrams, 
published by the Supernova Cosmology Project (hereafter 


SCP; Perlmutter et al. 


1997 


1998, [1999D and t he High- 


Redshift Supernova Search Team (hereafter HZT; Schmidt 


et al. 


1998; Riess et al 


1998). Data from both teams convinc¬ 


ingly reject the deceleration expected from an Einstein-de 
Sitter (EdS) Universe, and provide new evidence for a cos¬ 
mic acceleration in a low mass-density Universe consistent 
with a non-zero vacuum energy density. Evidence for spa- 
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199S|). 


Such important results demand excellent supporting ev¬ 
idence. In particular, it is necessary to examine the homo¬ 
geneity, environmental trends and evolutionary behaviour 
of SNela. It must be remembered that it is not yet clear 
whether these events are produced by a single mechanism 
in a unique physical environment. Systematic differences in 
the peak magnitudes between high and low-redshift SNe, 
arising either from evolutionary effects in the SNe progeni¬ 
tors or via subtle differences in the environments of the low 
and high-redshift galaxies in which SNe are produced, could 
mimic the cosmic acceleration without necessarily destroy¬ 
ing the small dispersion seen in existing Hubble diagrams. 

Differences between low and high-redshift SNela might 


199£; 

Hoflich et al. 200C; Dominguez et al, |200lD, different 

progenitor ages (Ruiz 

-Lapuente, Canal, & Burkert 

1997), 

greater amounts of dust in high-redshift environments, ei- 

ther in the host galaxy (Hatano, Branch, & Deaton 

1998; 

Totani & Kobayashi 

199t; Rowan-Robinson 2002|) 

or in 


the intergalactic medium (Aguirre 


1999 


Aguirre & Haiman 


200C), or a dependence of the SNe properties on host galaxy 


environments. Tests for analy sing SNe by the t - ype of their 


host galaxy were discussed by Perlmutter et al.| 


after P97), and continued by Perlmutter et al.| (199 



after P99). This initial classification based on host galaxy 
spectra revealed no changes in the properties of SNe located 
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in E/SO and spiral hosts, though only 17 hosts could be clas¬ 
sified and little morphological information was available. 

Continning these studies, via an examination of the 
dependence of SN properties on host galaxy morphology 
via Hubble Space Telescope (HST) imaging and Keck spec¬ 
troscopy, forms the basis of the present analysis. Although 
SNela can occur in all types of galaxies, disc and spheroidal 
stellar populations will represent different star-formation 
histories, metallicities and dust content. Thus we might ex¬ 
pect that SNela progenitor composition and peak magni¬ 
tudes or light-curve properties could be affected accordingly. 

Much of the previous work examining the environments 
of SNe la, ha,s been conducted using low-redshift sa.mnles. 


rehnement suggests that the earlier |Riess et al.| ( |l99S| ) trend 
may partly be due to the mix of E/SO and spiral galaxies, 
with (intrinsically fainter) SNe hosted by E/SO galaxies be¬ 
ing discovered at larger galactic distances when compared 
to SNe located in spiral galaxies. No significant trend is 
seen in SN parameters with galactocentric distance for the 
E/SO subset of the SN host population. As elliptical galax¬ 
ies show radial gradients in their stellar properties arising 
mostly from trends in metallicity rather than age, it is con¬ 
cluded that the intrinsic diversity in SN parameters more 
likely arises from stellar pop ulation differences r ather than 
progenitor metallicity alone (Ivanov et al. 2000). T he only 


relevant study for high- redshift samples is that of Howell 


with or ly limited work at high-redshift. However, in these _ Wang, fc Wheeler (200C), based on SCP lAU circulars and 


low-redshift studies, two interesting trends have already 
been claimed, both of which are relevant for this present pa¬ 
per. The first correlates the f orm of the SN Uqht-cur v e with 
the host-galaxy properties. Filippenko & Sargent (198£) 
originally found that the initial post-maximum decline rate 
in the light-curve in the B-band may have a smaller disper¬ 
sion among SNe located in elliptical galaxies. Usin g SNela 
disc overed v i a the Calan-Tololo (C -T) SN survey (Hamuy 
et al, 1996a), Hamuy et al. (1996b) also note that the post¬ 
maximum decline rate correlates with morphology, in the 
sense that SNe with more rapidly declining light-curves (i.e. 
intrinsically dimmer SNe) occur in earlier-type galaxies. In 


an ind ependent study of the CfA SN sample, Riess et al 


(1999) find a similar result. The various low-reds hift SN 
sam ples wer e later combined in a single analysis by Hamuy 
et al.^P 000| ). They conhrm the hndings of the C-T and CfA 
studies, and also note that intrinsically bright er SNe tend 
to occu r in bluer stellar e nvironments (see also Branch, Ro- 
manish n, & Baron 1996), and suggest that the intrinsically 


brightest SNe occur in the least l uminous galax ies. These 
effects were further investigated by Howelj (^001 ), who Hnd 
that in the local Universe over-luminous SNe arise in spiral 
and other late-type systems, whilst spectroscopically under- 
luminous SNe tend to occur in older (E/SO) systems. 

The second suggested correlation relates the photomet¬ 
ric pro perties of a SN to its physi c al loc ation in the host 
galaxy. Wang, Hoeflich, & Wheeler (1997), using SNe from 
the C-T survey as well as 11 other well-studied nearby SNe, 
found that SNe properties apparently vary with projected 
distance from the centre of the host galaxy. They find that 
SNe located more than 7.5 kpc from the centre show around 
3-4 times less scatter in maximum brightness than those 
within that projected radius, with E/SO galaxies dominat¬ 
ing the sample at high separations. Monte-Carlo simulations 
exploring the effects of dust on ra dial distributions of SNela 
( Hatano, Branch, fc Deaton 1998 ) similarly predict that SNe 
at large projected radial distances should be practically un¬ 
extinguished, whilst those at small projected distances could 
suffer significant amounts of extinction. The conclusion is 
that SNe at higher projected distance are a more homoge¬ 
neous group and should be better distance indicators. 


In a combined C-T/CfA sample, Riess et al. (199£) 
Hnd that SNela with faster decline rates (i.e. intrinsically 
fainter SNe) occur at larger projected dis tances from the 
host galaxy ce ntre. The Hamuy et alj (bOPCj) sa mple has also 
been used by Ivanov, Hamuy, fc Pinto (200C) for a similar 
purpose, though importantly using de-projected galactocen¬ 
tric distances rather than simple projected distances. This 


finding charts. These authors examine the radial distribution 
of the SCP SNe from their host galaxies, and find that they 
are distributed similarly to those from low-redshift CCD- 
based samples. 

The major difficulty in interpreting these trends lies 
in disentangling all of the numerous variables (morphol¬ 
ogy, radial position, luminosity, metallicity, dust content and 
star-formation history) which define the stellar populations 
of galaxies. Though no complete physical picture has yet 
emerged as to how the underlying stellar population might 
influence SN properties in t he manner propose d, one possi¬ 
ble hypothesis is as follows (Umeda et al. 1999). 

SNela most likely originate via the thermonuclear ex¬ 
plosion of a carbon-oxygen (CO) white dwarf (WD) in a 
binary system, with the light-curve powered by the decay 
of ®®Ni. Hence, WDs with a larger ^^C mass fraction will 
generate more ®®Ni and consequently brighter SNe. As the 
oldest progenitors are likely to have the lowest mass com¬ 
panions, less mass can be transferred to the primary WD, 
and the resultant reduced final ^^C abundance leads to a 
smaller ®®Ni mass, and therefore fainter events. In this sim¬ 
ple picture, one expects firstly to see more variation in SN 
properties in disc-based systems (where the age of the stellar 
population is more diverse) than in spheroidal galaxies, and 
secondly to see fewer over-luminous events in (presumably) 
older E/SO systems. 

So far as the implications for the use of SNe la as cosmo¬ 
logical probes are concerned, we are, of course, interested in 
possible trends in the above correlations with redshift. Our 
aim is to include SNe located at high redshift which are used 
in existing Hubble diagrams to constrain the cosmological 
parameters. Environmental effects on the determination of 
the cosmological parameters can then be more easily ex¬ 
plored. An outline of the paper follows. In §2 we discuss the 
low and high-redshift SN samples used in our analysis and 
describe the method utilised to align them on to the same 
calibration system. §3 introduces our new dataset, based 
primarily on HST imaging and Keck spectroscopy, used to 
investigate the properties of the host galaxies of the high- 
redshift sample. We discuss the classification system for the 
hosts in §4, and present our analyses in §5 and discussion in 
§6. We present our conclusions in §7. 


2 THE SUPERNOVA SAMPLES 

In order for us to explore the cosmological consequences of 
possible correlations between the SN properties and those 
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of the host galaxies, we need samples of SNela spanning the 
full range in redshift in the Hubble diagram analyses. 

For the distant SNela, we perform the analysis on the 
published sample of 42 SNe discovered by the Supernova 
Cosmology Project (SCP; see P99). All SNe were discov¬ 
ered by the SCP and details of the various procedures and 
analysis techniques, as well as spectroscopic verification, can 
be found in P97 and P99. 

The bulk of our sample overlaps with the P99 set but 
this overlap is not complete by virtue of the random man¬ 
ner in which HST snapshot observations are selected. For 
example, only 30 of the 42 have HST imaging providing 
host galaxy morphologies. P99 further restrict this sample 
of 42 SNe to one of 38 high-redshift objects in the primary 
cosmological fit (‘fit-C’ in P99) by excluding those SNe con¬ 
sidered ‘likely reddened’ or which are ‘outliers’ in terms of 
their stretch of the light-curve time-scale values. Where ap¬ 
propriate in this current analysis we also note the effect of 
excluding these 4 objects. 


At lower redshifts, we use the Hamuy et al. (1996a) SNe tude (rrTg'“) via 


sample and add selected SNe from the [Riess et^ 


199f) 


local sample. Our principal criteria in selecting SNe from 
these local samples remains as in P99, namely (i) each SN 
must have been observed prior to 5 days after maximum 
light so as to reduce any extrapolations necessary to fit peak 
luminosities, and (ii) that each SN has cz > 3000 km s“^ to 
control the impact of peculiar velocities. These requirements 
mean that 18 of the 29 C-T SNe (see P99 for details) and 
12 of the 22 Riess et al. (1999) CfA SNe^ are included in 
our low-redshift sample of 30 objects. As in P99, the two 
large stretch/residual outliers in the C-T sample (1992bo 
and 1992br, both located in E/SO systems) are excluded in 
the analysis, as well as two from the CfA sample (1994t 
and 199^e) _ One further SNe whose light-ciir\rp fit results in 


a nnrealsona.hly high stretch value is also excluded (199fia,h 
from the CfA sample). This leaves 25 low-redshift SNe in 
our primary analysis sample (16 from the C-T survey, and 
9 from CfA). 

It is important to ensure that the corrected peak mag¬ 
nitudes of the SNe drawn from the different samples are 
derived in a consistent manner. Raw observed SN pe ak mag- 


nitudes 


hav e an r.m.e. disp e rsion of < 0.30 mag ( e .g. Hamuy 


by using an empirical relationship between the light-curve 
shape and the raw observed (intrinsic) peak magnitude. Dif¬ 
ferent methods have been used to perform this correction 
and we are now concerned with selecting a single methodol- 
ogy for our c o mbin ed sample. 

Phillips (1993) showed that the absolute magnitudes 


of SNela are tightly correlated with the initial decline rate 
of the light-curve in the rest-frame B-band. They used a 
linear relationship to correct the observed SN peak magni¬ 
tudes, parameterized via Amis{B), the decline in B-band 
magnitudes over the first 15 days after maximum light - 
intrinsically brighter SNe have wider, more slowly d eclin¬ 


ing lig ht-curves. Further refinements by Hamuy et al. (1995 


1996b) can reduce the observed scatter to < 0.17 mag, and to 


1 SN1994m, SN1994s, SN1994t, SN1995e SN1995ac, SN1995ak, 
SN1995bd, SN1996ab, SN1996bv, SN1996bl, SN1996bo, SN1996c 


0.11 mag i f reddening correcti ons based on the SNe colours 
are made (Phillips et al 1999). 


An alternative technique, t he Multi-colo r Ligh t Curv e 


Shape (MFCS), introduced by Riess et al, (1995, 1996), 
adds or subtracts a scaled ‘correction template’ to a stan¬ 
dard light-curve template which creates a family of broader 
and narrower light-curves. A simple linear (or non-linear) 
relationship between the amount of the correction template 
added and the absolute magnitude of the SN can be used to 
correct the SN magnitude, and results in a similar smaller 
dispersion in the final peak magnitudes, particularly if multi¬ 
colour light-curves are utilised in the htting procedure. 

A third approach, and the one used in this paper, is 
to parameterize the SN light-curve time-scale via a simple 
stretch factor, s, which linearly stretches or contracts the 
time axis of a template SN light curve around the time of 
maximum light to best-fit the observed light curves of every 
SN (Perlmutter et al. 1996, 97, 99). A SN’s ‘corrected’ peak 
magnitude (mg””) is then related to its ‘raw’ peak magni- 


mg*”” = m”a"'" -b a(s - 1) 


( 1 ) 


where s is the stretch of the light-curve time-scale and a re- 
lat es s to the size of the m agnitude correction. As discussed 
in Goldhaber et al. (2001), the stretch-correction technique 
is at least as good at fitting existing data as any other single 
parameter method can be. It is important to note that, as 
we will show in later sections, not applying any correction 
to the peak magnitudes makes no signihcant difference to 
our conclusions. 

Each SN (at both high and low redshift) was fitted to 
determine the stretch factor s, using the published light- 


curves of the low-redshift samples (Hamuy et al. 1996a; Riesi 


et al 


1999) and the ‘exponential’ SCP light-curve template 


used in P99. The new fits for the CfA SNe are very sim¬ 
ilar to those of the C-T SNe. Each SN peak magnitude is 
corrected to the value appropriate for a fiducial light-curve 
corresponding to that of the template SN with s = 1. Fi¬ 
nally, every stretch-corrected SN B-band peak magnitude is 
also corrected for galactic extinction us ing the colour excess 


at the Galacti c coordinates of each SN (schlegel, Finkbeiner 


et al. 11961]), but this dispersion can be significantly reduced _ fc Davis 1^98). 


3 NEW DATA 

In this section we introduce the new dataset we have assem¬ 
bled for the host galaxies of previously-discovered SCP high 
redshift SNela. Comparable data on the lower redshift SNe 
are taken from the articles referenced in Section ^ 

The new data have been gathered through three ob¬ 
servational programs: (i) HST observations which provide 
spatially-resolved imaging data for morphologies, (ii) high- 
resolution ground-based spectroscopy which characterises 
the star-formation properties, and (iii) a variety of ground- 
based imaging and low-resolution spectroscopic data which 
serves to add to our capabilities of classifying the host galax¬ 
ies. A component of the latter survey includes the original 
SCP ground-based photometry of the host galaxies drawn 
from the reference frames (Aldering et ah, in prep.). 
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3.1 HST Imaging 


Our HST imaging data comes from two sources. The pri¬ 
mary source is a dedicated Space Telescope Imaging Spec¬ 
trograph (STIS) snapshot program (ID: 8313, 9131; PI: El¬ 
lis) of the host galaxies obtained during cycles 8 and 10. The 
second source of HST data arises via several (some serendip¬ 
itous) Wide Field Planetary Camera 2 (WFPC-2) observa¬ 
tions of a small number of the galaxies. A complete listii^ 
of the host galaxies observed with HST is given in Table hi 
Examples of the STIS images can be found in Fig. 


In the case of the STIS data, each targeted galaxy was 
imaged in the 50CCD (clear aperture) mode, which approx¬ 
imates a broad V -\- R-\-1 band-pass, in a 3 x 434 s dithered 
series. Basic image processing (bias subtraction, dark sub¬ 
traction and fiat-fielding) was performed using the stan¬ 
dard STIS reduction pipeline from within iraf. Cosmic ray 
removal and final image construction was then performed 
using a ‘variable pixel linear-reconstructi on’ technique us¬ 
ing t he DITHER-II and drizzle packages (Fruchter & Hook 


1997). The output pixel size from this process was 0.025 


arcsec across, and a value pixfrac = 0.6 was used. 


As well as this dedicated STIS snapshot program, we 
use WFPC-2 data for one host from another SCP HST pro¬ 
gram (ID:7590, PI: Perlmutter), and also searched the HST 
archive for additional HST observations of the hosts, finding 
three cases independently observed by HST /WFPC-2 (see 
Table and references in P97 and P99). The WFPC-2 data 
sets were again combined using the dither-ii package used 
to reduce the STIS data. 


The location of the faded SNe cannot be predicted with 
certainty on the HST (in particular STIS) images due to 
the pointing uncertainties generic in snapshot observations. 
In some cases where more than one galaxy is located near 
the centre of the HST pointing, this could confuse the cor¬ 
rect assignment of the host galaxy. To accurately locate the 
position of the SNe on the HST (STIS or WFPC-2) images 
therefore requires careful comparisons with existing ground- 
based data taken in good seeing conditions. We performed 
a cross-correlation between the HST images and the exist¬ 
ing calibrated ground-based images (cf. Aldering et ah, in 
prep.). In these images, the SNe are visible, and hence the 
exact position can be determined in relation to surround¬ 
ing objects. For the STIS field-of-view of approximately 
52" X 52", typically 6 — 8 compact stellar-like objects were 
common to both the deep ground-based and the STIS data. 
We calculate the relationship between the HST and ground- 
based pixel coordinate systems using GEOMAP and related 
tasks within iraf, and are then able to transform the pixel 
position of the SNe from the ground-based coordinate sys¬ 
tems to the HST frame of reference. 


We estimate the uncertainties in the derived positions 
using a simple ‘bootstrap’ technique. We calculate the loca¬ 
tion of the SN using all the common objects, and then per¬ 
form the same calculation but omitting one common object 
each time. The resultant variation in the calculated location 
of the SN gives approximate uncertainties in the derived lo¬ 
cation using all the common objects. The process typically 
provides SN positions which are accurate to ~ 0.15-0.4 arc- 
sec; the uncertainty is dependent on the number of objects 
in common between the two images and the seeing of the 


ground-based data. This accuracy is, in all cases, sufficiently 
precise to unambiguously identify the host galaxy. 

This identification technique also enables us to derive 
projected radial distances from the SNe locations to the cen¬ 
tre of their host galaxies for a sub-sample of the SNe. The use 
of spatially resolved HST images has the advantage of be¬ 
ing able to precisely locate the centre of each galaxy reliably, 
something which is more difficult when using ground-based 
images. The uncertainty therefore derives entirely from the 
accuracy to which we can locate the SN position. In 20 of the 
36 SNe with morphological information from HST images, 
we are able to calculate the SN location to an accuracy of 
< 0.3 arcsec; these data will be discussed further in a later 
article. 

We convert these distances in arcseconds to projected 
distances in kpc using the angular diameter-distance re¬ 
lation; where necessary in these calculations we assume a 
Hm = 0.28, Ha = 0.72 and Hq = 70kms“^Mpc“^ cosmol¬ 
ogy- 


3.2 Keck spectroscopy 


The pre-existing ‘discovery’ host galaxy spectra were taken 
during the period when the SN light dominated that of the 
host galaxy. The principal goal of those spectroscopic cam¬ 
paigns was to determine a precise type and redshift for the 
SN candidate. As galaxies typically have narrow emission 
and absorption line features, their presence can often be de¬ 
tected even at maximum SN light; however, the contamina¬ 
tion by the SN usually makes any detailed study of the host 
galaxy properties very difficult. We have therefore performed 
a second round of optical spectroscopy for a sub-sample of 
the host galaxies imaged by HST using the Echellette Spec¬ 
trograph and Imager (ESI) on t he Keck-II 10 m tel escope 
(details on ESI can be found in sheinis et al 2002). The 
goal of this component is to considerably improve the diag¬ 
nostic spectroscopy for the host galaxies. 

ESI can operate in three observing modes (imaging, low 
dispersion and echellette); each galaxy was observed in the 
echellette and imaging configurations. The echellette mode 
has a continuous wavelength coverage of ~ 4000-11000 A, 
spread across ten orders with a constant 11.4 kms“^ pixel” ^ 
dispersion. All spectroscopic observations were performed 
using a slit of width 1 arcsec. The position angle was calcu¬ 
lated from the HST images to ensure the slit passes through 
both the location of the (now faded) SN and the host galaxy 
centre. Table contains a log of the observations. 

The primary data reduction steps were performed us¬ 
ing the MAKEE data reduction software package written by 
T. Barlow^. This package produces a bias subtracted, flat- 
fielded, wavelength calibrated and sky subtracted, but not 
order-combined, spectrum of each object by using a master 
trace of a bright star and shifting this trace to the position 
of the object to be extracted. A master bias frame and in¬ 
ternal (quartz lamp) flat-field were created on every night 
of observation and used for that nights data; the wavelength 
calibrations were performed using HgNeXe and CuAr lamps 
bracketing each science observation. 


^ MAKEF, is availa.hle from 

http: //spider.ipac.caltech.edu/staff/tab/makee/index.html 
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Table 1. HST observing log for the host galaxies of the P99 high-redshift SNe discovered via the SCP collaboration. 


SN Name 

Program ID 

Configuration 

Exposure time (s) 

Type“ 

Notes 

1992bi 

8313 

STIS/50CCD 

3 X 434 

2 

Sm/Irr 

1994al 

8313 

STIS/50CCD 

3 X 434 

0 

SO+comp 

1994h 

7786 

WFPC-2/606W 

1 X 100,1 X 40 

1 

s/comp (see P97) 

1994an 

8313 

STIS/50CCD 

3 X 434 

2 

Scd 

1994am 

5378 

WFPC-2/814W 

3 X 2100 

0 

E (see P97) 

(1994am 

9131 

STIS/50CCD 

3 X 434 

0 

E) 

1994g 

8313 

STIS/50CCD 

3 X 434 

2 

Pec/Irr 

1994f 

9131 

STIS/50CCD 

- 


Guide star acquisition failed 

1995aq 

9131 

STIS/50CCD 

- 

- 

No observations 

1995ar 

8313 

STIS/50CCD 

3 X 434 

1 

Sa 

1995az 

8313 

STIS/50CCD 

3 X 434 

1 

Sc 

1995ay 

9131 

STIS/50CCD 

- 

- 

No observations 

1995at 

8313 

STIS/50CCD 

3 X 434 

1 

Sa 

1995ba 

9131 

STIS/50CCD 

3 X 434 

0 

E/SO 

1995ax 

8313 

STIS/50CCD 

3 X 434 

0 

E/SO 

1995aw 

8313 

STIS/50CCD 

3 X 434 


Low S/N 

1995as 

8313 

STIS/50CCD 

3 X 434 

2 

S/Irr 

1996cg 

8313 

STIS/50CCD 

3 X 434 

1 

Sab 

1996cf 

8313 

STIS/50CCD 

3 X 434 

2 

Scd 

1996ci 

9131 

STIS/50CCD 

3 X 434 

2 

Scd 

1996cm 

9131 

STIS/50CCD 

3 X 434 


Low S/N 

1996cl 

5987 

WFPC-2/814W 

6 X 2600 

0 

SO (see P99) 

(1996C1)'’ 

(7372) 

(WFPC-2/814W) 

(6 X 1100) 

(0) 

(SO) 

1996ck 

8313 

STIS/50CCD 

3 X 434 

0 

E/SO 

1996cn 

8313 

STIS/50CCD 

3 X 434 

0 

SO 

1997i 

8313 

STIS/50CCD 

3 X 434 

1 

Sbc 

1997h 

8313 

STIS/50CCD 

3 X 434 

0 

E 

1997g‘= 

9131 

STIS/50CCD 

3 X 434 

2 

Merging pair 

1997f 

8313 

STIS/50CCD 

3 X 434 

1 

SBb 

1997j 

9131 

STIS/50CCD 

1 X 434 

0 

E/S0;lnterrupted exposure 

19971 

8313 

STIS/50CCD 

3 X 434 

2 

Irr 

1997o 

9131 

STIS/50CCD 

- 

- 

Guide star acquisition failed 

1997p 

9131 

STIS/50CCD 

- 

- 

No observations 

1997s 

9131 

STIS/50CCD 

- 


No observations 

1997n 

9131 

STIS/50CCD 

3 X 434 

2 

Irr; poss. merger 

1997q 

9131 

STIS/50CCD 

3 X 434 

0 

so 

1997k 

8313 

STIS/50CCD 

3 X 434 

2 

Irr 

1997r 

9131 

STIS/50CCD 

3 X 434 

1 

Sab 

1997ac 

8313 

STIS/50CCD 

3 X 434 

0 

E/comp;low S/N 

1997aj 

8313 

STIS/50CCD 

3 X 434 

1 

Spec 

1997ai 

9131 

STIS/50CCD 

3 X 434 


Low S/N 

1997af 

8313 

STIS/50CCD 

3 X 434 

2 

Irr 

1997am 

9131 

STIS/50CCD 

3 X 434 

2 

Soec/Irr 

1997ap 

7590 

WFPC-2/814W 

3 X 700 

- 

Host not visible (see 

Perlmutter et al, 

1998) 


Notes: 

“ where: ‘0’ indicates E/SO type, T’ early-type spirals, and ‘2’ late-type spirals/irregulars 
** 1996cl has been observed twice by HST /WFPC-2, prior to and during the SN event (see e.g. 

1997g was also observed as part of 8313; however, this observation failed during guide star acquisition. 


Gal-Yam, Maoz, fc Sharon| 2002 


On occasion, the automatic profile locating routine in et al, 1987; Krisciunas 199C). Though an absolute flux cali- 


MAKEE failed for some orders or for those host galaxies which 
were spatially resolved. In these cases, the object profile lo¬ 
cation on the slit and the profile width were determined by 
hand and the spectrum extracted with these parameters. 

In order to combine the ten orders into a single con¬ 
tinuous spectrum, we performed flux calibration of the MA- 
KEE-extracted spectra using standard iraf routines using a 
flux standard observed on the night of observation at a sim¬ 
ilar airmass to the scie nce target; su i table stan dards were 
drawn fro m the lists of Massey et al. (1988) and Massey & 
Gronw^ (1990). We also correct for the extincti on proper- 
ties of the Mauna Kea summit atmosphere (e.g. Krisciunat 


bration is not possible with a slit of finite width, a relative 
flux calibration between spectral regions is more reliable. 
The primary uncertainty here arises from variations in the 
seeing as a function of wavelength, where the FWHM of an 
image varies as We estimate the effect to be at the 

^ 5 per cent level between Ha and H/3, the primary emission 
lines that interest us in this analysis. 

After flux calibration, each order was examined and the 
wavelength range offering the best signal-to-noise (S/N) de¬ 
termined. These wavelength ranges were then median com¬ 
bined to form a final flux-calibrated spectrum from each in¬ 
dividual exposure, and finally the individual exposures me- 
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1995axE/S0 1996ck E/SO 



1997aj Sa/pec 1996cg Sab 





, 'Q^; .: 

1997h E/SO 1994al SO 



1997fSBb 1992bi Sm/Irr 




1995azSc 1997i Sb/c 


1994an Scd 1995as S/irr 


Figure 1. A selection of the SNela host galaxies imaged using STIS on the HST, demonstrating the range of galaxy types and the 
quality of data produced by the program. In each image, the circle indicates the location of the (faded) SN derived via comparisons with 
existing ground-based data. 


Table 2. The ground-based observing log for the new data on the host galaxies of the SCP high-redshift SNe. 


SN Name 

Spectroscopic 
observation date 

Spectroscopic 
exposure time (s) 

Imaging 

observation date 

Imaging data 

source 

1992bi 

2001 May 21 

2 X 1200 

2001 May 21 

ESI 

1994g 

2000 Feb 10 

3 X 1200 

2000 Feb 10 

ESI 

1994al 

2000 Nov 24 

2 X 1800 

2000 Oct 31 

COSMIC 

1994an 

2000 Nov 24 

2 X 1800 



1995ar 

2001 Nov 13 

2 X 1800 



1995az 

2001 Oct 20 

2 X 1500 

2000 Oct 31 

COSMIC 

1995ax 

2000 Nov 22 

2 X 1800 

2000 Oct 31 

COSMIC 

1996cg 

2000 Feb 10 

3 X 1200 

2000 Feb 11 

ESI 

1997i 

2000 Feb 10 

2 X 1200 

2000 Feb 10 

ESI 

1997h 

2001 Oct 23 

2 X 1500 



1997f 

2000 Feb 11 

3 X 1200 

2000 Feb 10 

ESI 

1997aj 

2001 May 22 

2 X 1200 

2001 May 22 

ESI 

1997g 



2000 Feb 10 

ESI 


dian combined to form a final output spectrum for each ob¬ 
ject. 

The spectra were examined using the SPLOT facility in 
IRAF. For each output spectrum a redshift was measured 
via visual inspection, typically using Call H,K absorption or 
[O ll]/Balmer emission-line features. We find excellent agree¬ 
ment with the earlier redshifts reported in P97/P99 in all 
cases. Fig. ^ demonstrates a range of the host galaxy spec¬ 
tra observed; comments on each spectrum can be found in 
Table ^ For the emission line spectra, we also attempt to 
measure the strength of the Ha and H/3 emission lines, where 
they are present and are not contaminated by poor residual 
sky subtraction from strong sky lines. These measures were 
also made usliig sPLU't. - 


in a host galaxy, as the reddening nature of the dust makes 
the observed ratio of the fluxes of two emission lines (Robs) 
differ from the intrinsic ratio (i?int). The colour excess of the 
ionised gas, E{B — F)g, can then be estimated via 

log (//obs ///int) 


E{B-V)^ = 


( 2 ) 


0.4[fc(AH/3) — fc(AHa)] 

where fc(A) is the form of the total-to-selective extinction 
of the gas for a particular extinction law, and Ah /3 and 
Ana are the rest-frame emission wavelengths of H/3 and Ha 
respectively. Ass uming case-B rec ombination, R^nt = 2.87 
for Ha/H/3 (e.g. Osterbrock 1989). The reddening parame¬ 


ter Ay can be calculated via Ay ~ E{B — V)„R v, where 
Ry ~ 3.1 for a Milky Way-type extinction law (e.g. |Cardelli 


et al. 


198C 


Pei 


1992). Where available, these Ay parame- 
The Ha and H/3 lines give an indication of the extinction ters are also listed in Table H also listed are values where the 
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J-l 

CD 



SN 1997h 
z = 0.526 
HST type: E 
Spectral type: 0 



SN 1994al 
z = 0.421 
HST type: SO 
Spectral type: 0 



SN 1995ar 
z = 0.465 
HST type: Sa 
Spectral type: 1 


X 

o 

o 



SN 1994an 
z = 0.378 
HST type: Scd 
Spectral type: 2 


- 


' - 




on 

H /3 OIII 


Ha 

SK 

Y : 

: 





Am 


5290 5300 5310 5320 6900 6950 7000 7050 7100 71EBB00 9350 9400 


Wavelength (A) 


SN 1994g 
z = 0.425 

HST type: Pec/Irr 
Spectral type: 2 


Figure 2. Example spectra taken from the ESI spectroscopic campaign demonstrating the range of spectral types and the quality of the 
ESI data obtained. 


extinction due to the Milky-Way has been removed. These 
Av parameters have been corrected for stellar absorption 
affecting the Ho? and H/3 lines assuming stellar absorption 
equivalent widths of 2A in absorption on each. We make no 
corrections for the seeing variations between the two wave¬ 
lengths of interest, but note that as seeing typically improves 
with wavelength (and hence slit losses will be smaller), our 
measured Hq to H/3 ratios will effectively be upper limits on 
the intrinsic ratios. 

One further subtlety is that these values of Ay are mea¬ 
sured for the extinction on the nebular emission lines origi¬ 
nating in star-forming H ll regions. This may not, and indeed 
probably will not, be representative for regions of the galaxy 
containing older stellar populations where SNIa may orig¬ 
inate; indeed, it has been shown that the continuum emis¬ 
sion from stars in starburst galaxies is often less obscured 
than the line emissio n from the gas (e.g. ^hnelhet^ 


1988 


Calzetti et al. 1994; Mas-Hesse & Kunth ~ p.999D . An esti¬ 
mate o f this extinctio n can be derived using the prescrip¬ 
tion of Calzetti et al ( 200C ), who relate the colour excess 
of the stellar continuum, E{B — V)b, to E{B — I/)g using 
E{B - V)s = 0A4:E{B - V)g. 


We discuss the utility of the spectroscopic data further 
in Sections ^ and 


3.3 Ground-based Imaging 

Many of the original SN search images detect the host 
galaxy, and photometric measures of these hosts are de¬ 
scribed and analysed in a companion paper which includes 
the reference and late-time images free from SN contam¬ 
ination (Aldering et ah, in prep). As ESI offers an direct 
imaging capability, we took the opportunity to gather addi¬ 
tional photometric measures when undertaking the spectro¬ 
scopic component described above. Additional photometry 
was obtained with the COSMIC instrument at the Hale 200- 
in telescope at the Palomar Observatory. An observing log 
of all the new imaging observations obtained in addition to 
that described in Aldering et al. is given in Table 

The galaxies were observed in filters which most closely 
matched those filters used in the original SNe search pro¬ 
gram and in Aldering et al., the standard Cousins R and I 
band filters (corresponding to rest-frame bands B and V for 














































The Hubble Diagram of Type la Supernovae as a Function of Host Galaxy Morphology 9 


most of our sample). In the case of ESI we used standard 
R and Gunn i filters ( Thuan fc Gunn| 1976), and for GOS- 
MIG we used Gunn r and i. The new imaging data were re¬ 
duced using standard techniques: overscan subtraction, bias 
subtraction and fiat-fielding were performed using calibra¬ 
tion data taken on the night of observation. The target host 
galaxy was identified using the original SN detection images, 
and raw magnitud es measured using the SEXTRACTOR pro¬ 
gram version 2.21 (Bertin & Arnoutt 1996). We used the au¬ 
tomatic aperture magnitudes determined by SEXTRACTOR, 
whic h are similar to Kron’s ‘first moment’ algorithm (Kron 


198C). As the R and I images were taken consecutively on 


each night of observation, variations in the seeing between 
the R and I band data are negligible. 

All the imaging dat a were taken on photometric n ights 
and standard star fields ( |!^andolt 1992; Jorgensen 1994) were 
observed at various airmasses to derive the GCD zeropoints 
and variation of extinction with airmass in order to calibrate 
the observed raw magnitudes on to the a-Lyrae system. Fi¬ 
nally, each integrated magnitude was correc ted for galactic 


extinc tion using the colour excess maps of Schlegel et al 
( [1998| ). 

Seven of the 11 host galaxies imaged with ESI i?-band 
filter overlapped with the sample of Aldering et ah, en¬ 
abling us to determine the photometric precision of the joint 
dataset. We directly compare the i7-band magnitudes de¬ 
rived from the two datasets, and find an excellent agree¬ 
ment; the mean absolute difference between the two datasets 
is only 0.06 mag (in the sense that the ESI magnitudes are 
fainter), with the largest discrepancy 0.12 mag. This discrep¬ 
ancy arises principally from the different magnitudes mea¬ 
sured: Aldering et al. define their magnitudes to include all 
the light from within two Petrosian radii of a galaxy, whilst 
the definition of the magnitudes in SEXTRACTOR leads to a 
slightly smaller aperture. These differing definitions will not 
affect the calculated colours unless the galaxies under study 
have marked colour gradients in their outermost regions. In 
the future analyses in this paper where we use these colours 
to help classify the hosts galaxies, we use data in the pref¬ 
erence: i) Keck-ESI R,i, ii) Gousins R,I (Aldering et al, in 
prep.), and iii) COSMIC r,i. 


4 HOST GALAXY CLASSIFICATION 

We now turn to the question of characterising each host 
galaxy according to the data presented above. In the case of 
the high redshift SNe, we classified each host galaxy accord¬ 
ing to one of three broad categories: spheroidal (E/SO, type 
0), early-type spiral (Sa through Sbc, type 1) and late-type 
spiral (Scd through to irregular types, type 2), based on a 
combination of the three diagnostics available to us. 

In order of preference in the classification scheme, these 

are: 

• HST morphology. The morphologies of the galaxies 
were visually classified by two of the authors (RSE and MS). 
The morphology of each HST-observed host is given in Ta¬ 
ble examples of each classification can be found in Fig. 

• Ground-based spectrum. The ESI spectrum (or if not 
available a pre-existing SN discovery spectrum with suffi¬ 
cient S/N) was classified according to the presence or ab¬ 
sence of nebular emission lines and the prominence of the 


discontinuity around Gall H and K. Our scheme runs from 
0 (no emission lines, strong break) to 2 (clear strong H ll- 
like spectrum, no break). Comments on each object can be 
found in Table ^ example spectra are shown in Fig. |^. 

• Ground-based inte grated colour. W e use a set of tem¬ 
plate SEDs (taken from Poggiant: 1997) to estimate the ap¬ 
proximate galaxy class based on the observed R-I colour and 
known redshift. The appropriate R (or r) and I (or i) filter 
response curves are convolved with each SED template at 
the host redshift, with the galaxy assigned a ‘colour-type’ 
based on the nearest matching galaxy SED. We do not ac¬ 
count for redshift evolution in the galaxy SEDs in this pro¬ 
cess. 


The question now arises as to how to combine these clas¬ 
sifications, recalling that not all diagnostics are available for 
each host galaxy. We assigned a confidence of either 0 (prob¬ 
able) or 1 (secure), based on the agreement and availability 
of the various measures. This leads to two different samples 
of data - a ‘securely’ classified sample (our primary analysis 
dataset) and a slightly enlarged, but less securely classified, 
sample. Note that even in the less secure sample, the dis¬ 
crimination between spheroidal and disc based systems is 
usually robust, even if the precise class in the case of spiral 
systems is not clear. The following scheme is used: 

• The HST imaging and ESI spectra; in all available cases 
these agree (Gonfidence ‘1’), 

• Clear HST morphology; no ESI spectrum (Confidence 

‘ 1 ’). 

• Ambiguous HST morphology; supporting evidence 
from pre-existing (SN discovery run) spectral type or colour 
(Confidence ‘1’), 

36 host galaxies (86 per cent of the total of 42 SNe 
in the P99 high-redshift sample) meet the above criteria. 
Three additional galaxies have been less securely classified 
according to: 

• No HST image; class taken from good quality ground- 
based imaging and spectroscopic data (Confidence ‘0’). 

• No HST image; class taken from good quality ground- 
based imaging or spectroscopic data (Confidence ‘0’). 

In total we therefore have 39 classified galaxies in the 
enlarged sample (93 per cent of the total P99 high-redshift 
objects). The classification statistics can be found in Table ^ 
and the agreement between the different assigned types for 
the hosts can be examined in Pig. ^ in the context of the 
adopted {R — I) colour-redshift relation for those galaxies 
for which R and 7-band measures are available. The clas¬ 
sifications are colour coded from red (E/SO) through green 
(early spiral) to blue (late spiral/irregular); black denotes 
those galaxies with no assigned class; these data can also be 
found in Table ^ 

As noted, the classifications of each object using the 
different techniques agree well; in particular, there are no 
red objects with emission line spectra or spiral morphol¬ 
ogy and no blue objects without emission lines which have 
a spheroidal HST-type. The greatest difficulty is in distin¬ 
guishing between type 1 and type 2 objects, where the R — 1 
colour is of only limited value, possibly due to the increased 
presence of dust or evolutionary effects in these spiral sys¬ 
tems when compared to E/SO objects; in these objects we 
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Table 3. ESI ground-based spectral classification 


SN Host 

Spectral type 

Ay (raw) ^ 

Ay (corrected)^ 

Gomments^ 

1992bi 

2 

0.11 

0.06 (0.03) 

Broad [Oil], H/3, [Oiii] 

1994g 

2 

0.64 

0.62 (0.27) 

s. [On], H/3, [Oiii] 

1994al 

0 

- 

- 

s. H,K, B4000 

1994an 

2 

1.08 

0.87 (0.38) 

[On], H/3, [Oiii], Ha 

1995ar 

1 

- 

- 

H,K, B4000, [On], H/3 

1995az 

1 

3.71 

3.16 (1.39) 

[On], H/3, Ha 

1995ax 

0 

- 

- 

Flat; w. H,K, B4000 

1996cg 

1 

- 

- 

w. [O ii] 

1997i 

1 

1.15 

1.01 (0.44) 

s. [On], H/3, [Oiii], Ha 

1997h 

0 

- 

- 

s. H,K, B4000 

1997f 

1 

- 

- 

w. [Oil], H/3, [Oiii] 

1997aj 

1 

- 

- 

[On], H/3, [Oiii] 


Notes: 

“B4000=4000A break, s.=strong, w.=weak 

^From Balmer decrement (corrected for stellar absorption), where available 

^Corrected for Galactic extinction, bracketed value is based on the Calzetti et al. (2000) prescription; see text. 



Redshift 


Figure 3. The colour-redshift relation for the subset of SN host galaxies with both R and I measures, together with the classification 
of the SN host galaxies from the three diagnostics available (the colour type, HST-type and spectral type). The colour shown is Galactic 
extinction corrected R-I; galaxies with colours measured based on Gunn filters are not sh own fthough sim ilar diagrams exist for these 
smaller subsets). The final adopted classification is also shown. Overlaid SEDs are from Poggianti (1997), ranging from elliptical (E) 
to starburst (SB) galaxies. In the classification process, no redshift evolution in the galaxy properties is assumed to occur. Internal 
extinction in the host galaxies is not accounted for on this plot. 
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Table 4. The classification of the P99 42 high-redshift SN host galaxies 


SN name" Redshift ESI ST^ Other ST^ Colour HST Final Classification 

Type Type Type security 


1992bi 


0.458 

2 

1 

2 

2 

2 

1 

1994al 


0.420 

0 

0 

- 

0 

0 

1 

1994h 

* 

0.374 

- 

0 

- 

1 

1 

1 

1994an 


0.378 

2 

2 

2 

2 

2 

1 

1994am 


0.372 

- 

0 

- 

0 

0 

1 

1994g 


0.425 

2 

2 

1 

2 

2 

1 

1994f 


0.354 

- 

1 

1 

- 

1 

1 

1995aq 


0.453 

- 

2 

2 

- 

2 

1 

1995ar 


0.465 

1 

- 

1 

1 

1 

1 

1995az 


0.450 

1 

1 

2 

1 

1 

1 

1995ay 


0.480 

- 

2 

2 

- 

2 

1 

1995at 


0.655 

- 

0 

1 

1 

1 

1 

1995ba 


0.388 

- 

1 

2 

0 

0 

1 

1995ax 


0.615 

0 

0 

1 

0 

0 

1 

1995aw 


0.400 

- 

- 

0 

- 

0 

0 

1995as 


0.498 

- 

0 

1 

2 

2 

1 

1996cg 

* 

0.490 

1 

- 

2 

1 

1 

1 

1996cf 


0.570 

- 

- 

2 

2 

2 

1 

1996ci 


0.495 

- 

- 

1 

2 

2 

1 

1996cm 


0.450 

- 

- 

2 

- 

2 

0 

1996cl 


0.828 

- 

0 

1 

0 

0 

1 

1996ck 


0.656 

- 

0 

1 

0 

0 

1 

1996cn 

* 

0.430 

- 

- 

1 

1 

1 

1 

1997i 


0.172 

1 

1 

2 

1 

1 

1 

1997h 


0.526 

0 

0 

1 

0 

0 

1 

1997g 


0.763 

- 

- 

2 

2 

2 

1 

19971 


0.580 

1 

1 

2 

1 

1 

1 

1997j 


0.619 

- 

0 

1 

0 

0 

1 

19971 


0.550 

- 

- 

2 

2 

2 

1 

19970 

* 

0.374 

- 

2 

2 

- 

2 

1 

1997p 


0.472 

- 

1 

1 

- 

1 

1 

1997s 


0.612 

- 

2 

- 

- 

- 

- 

1997n 


0.180 

- 

2 

2 

2 

2 

1 

1997q 


0.430 

- 

- 

1 

0 

0 

1 

1997k 


0.592 

- 

- 

2 

2 

2 

1 

1997r 


0.657 

- 

- 

1 

1 

1 

1 

1997ac 


0.320 

- 

- 

2 

0 

- 

- 

1997aj 


0.581 

1 

1 

2 

1 

1 

1 

1997ai 


0.450 

- 

- 

- 

- 

- 

- 

1997af 


0.579 

- 

2 

2 

2 

2 

1 

1997am 


0.416 

- 

2 

2 

2 

2 

1 

1997ap 


0.830 

- 

- 

2 

- 

2 

0 


Notes: 

"Starred SNe are excluded from the primary fit (‘fit-C’) of P99. 
^ST = Spectral Type. 


Table 5. Statistics on the classification of the P99 42 high-redshift plus the extended low-redshift sample SN host galaxies 

Host type Low redshift" High redshift" Total" 




Secure 

All 

Secure 

All 

E/SO 

3(5) 

9(9) 

10 (10) 

12 (14) 

13 (15) 

Early Spi 

14 (14) 

9 (12) 

9 (12) 

23 (26) 

23 (26) 

Late types 

8(8) 

14 (15) 

16 (17) 

22 (23) 

24 (25) 

Small Distances 

14 (14) 


19 (19) 


33 (33) 

Large Distances 

10 (12) 


7 ( 7) 


17 (19) 

All 

25 (27) 

32 (36) 

35 (39) 

57 (63) 

60 (66) 


11 


Notes: 

"Bracketed numbers include SNe excluded from our primary sample as stretch/residual outliers, or those likely reddened; see text and 
P99 for further details. 



12 M. Sullivan et al. 


mostly rely on the other two diagnostics. Thus it is impor¬ 
tant to recognise the most robust distinction in our survey 
lies between type 0 and the remainder. 

Classihcations for the low-redshift SN host galaxy sam¬ 
ple are taken from Hamuy et al. (2000). These authors give 
more than adequate morphological classes for all of the SNe 
in the low-redshift sample, and we simplify their system into 
the same three broader classes that we use for the higher 
redshift sample. Additional data on the projected distance 


( in ar csec) from the host galaxy is taken from Ivanov et al 


( 2000 ), which we convert into a projected distance in kpc. 

One potential concern with the consistency of the mor¬ 
phological typing might be a mismatch between the Hlters 
used for the imaging at low and high redshifts - a morpho¬ 
logical fe-correction. However, the breadth of the STIS hlter 
in use should minimize this effect. The low-redshift sample 
is ty pically typed using B or H-band images (see Hamuy 


et al, 2000), which at the mean redshift of our high-redshift 
sample would correspond to a i? or / passband, which falls 
within the V+B+I response of the STIS 50CCD; indeed this 
was one of the motivations for selecting the broad 50CCD 
observing mode. 


5 RESULTS 

With each galaxy classihed according to its morphological or 
spectral type, we now search for any environmental effects 
that may be present in the combined sample of high and low 
redshift SNe. This section presents our results in two broad 
categories. Firstly, we examine the Hubble diagram and its 
associated statistics as a function of the environment of the 
SN, both in terms of the host galaxy morphology and galac- 
tocentric distance. Secondly, we investigate possible correla¬ 
tions between the light curve parameters of the SNe and its 
environment similarly characterised. We discuss the broader 
implications of our findings in Section ^ 


5.1 Hubble diagrams as a function of SN 
environment 

We hrst examine the updated P99 Hubble diagram (hg. 2 
in P99), labelling each SN according to its host galaxy type. 
The upper panel of Fig. ^ shows this Hubble diagram in 
terms of the effective rest-frame B magnitude, i.e. corrected 
for the light-curve stretch-luminosity relationship of equa¬ 
tion (^), as a function of redshift for the 42 SCP high-redshift 
SNe, together with the 25 SNe from the low-z C-T/CfA sam¬ 
ples. A zoom to the high redshift sample is also shown. Those 
SNe excluded from the primary ht of P99 are marked. 

The SNe are colour-coded according to the host galaxy 
class. Less-securely classihed hosts are also marked. Those 
SNe with no classes are marked unknown and ignored in the 
subsequent analyses. The overlaid lines correspond to the 
apparent magnitude as a function of redshift for a standard 
s = 1 SNIa expected in an EdS Universe (Hm = 1, Ha = 
0) and a hat, A-dominated Universe with the best- htting 
cosmology of P99 (Hm = 0.28, Ha = 0.72; see P99, also jliess 


et al. 


1998). The lower panel in Fig. ^ shows the residuals 


from this best-htting line for each SN. 

We also perform additional htting procedures for the 
SNe residing in each galaxy class using a technique similar 


to that used in P99. We consider two cases - htting for 
both Hm and Ha independently, and htting for the case of 
a hat Universe (i.e. htting only for Hm or Ha and assuming 
Hm + Ha = 1). This involves solving the equation 

mg'''’ = -I- q(s - 1) = Mb + SlogDLiz', Hm, Ha), (3) 


where Bl = Hq di is the ‘Hubble-constant free’ luminosity 
distance, and Mb = Mb — 5 log Hq -|- 25 is the ‘Hubble- 
constant-free’ H-ban d absolute peak magnitude of a SNIa 
with s = 1 (see e.g. Uoobar & Perlmuttei 1995, and P99). 
This equation is solved for the parameters a, Mb, Hm and, 
for non-hat geometries, Ha, by least-squares i.e. minimizing 


(Hm, Ha, a, Mb) = 

(m^Br - [Mb + 5 log Dl jzi) - a{si - 1)])^ 

i—1 

As in P97 and P99, the small correlations between the high- 
redshift SNe (due to shared calibration data) are accounted 
for using a correlation matrix of the various uncertainties. 
We compare two numerical minimization techniques; a sim¬ 
ple Levenberg-Marquardt technique, and a more complex 
Gauss-Newton method, and hnd an excellent agreement be¬ 
tween the derived parameters in both cases. 

Details of the best-htting parameters are given in Ta¬ 
ble ^ The table also includes the ht obtained for all of the 
SNe with a host galaxy type (regardless of class) which pro¬ 
vides a valuable consistency check that our ‘typed’ sample 
is representative of the larger one studied by P99. The pa¬ 
rameters a and Mb are determined only for the P99 SNe 
sample and for the new enlarged sample, and the parame¬ 
ters from the latter ht adopted as hxed in the smaller sub¬ 
samples. (We study variations in Mb in detail in Section ^ 
whilst the a parameter has only a small impact at high- 
redshift due to the narrow stretch distribution seen, e.g. 
hg. 4 in P99). For these hts, the errors in the best-htting 
a and Mb parameters are carried through in the analy¬ 
sis. For a hat Hk = 0 Universe, Table ^ also lists the mean 
dispersions from the best-htting cosmologies for each sam¬ 
ple, as well as formal 1 -ct errors in Hm derived from the 
co-variance matrices of the various hts, assuming normally 
distributed errors. The htting procedure also generates two- 
dimensional conhdence regions, ■P(Hm,Ha) oc exp(— x^/ 2), 
shown in Fig. ^ for a selection of the SN sub-samples. The 
plot shows the 68% and 90% conhdence ellipses in the (Hm, 
Ha) plane together with the best-htting parameters for gen¬ 
eral and hat cosmologies. For those datasets where a and 
Mb are ht, as in P99 we calculate two-dimensional con¬ 
hdence regions by integrating over these parameters, i.e., 
P(Hm,Ha) = / / T’(Hm, Ha, AIb, a) dAIe da. This proba¬ 
bility space also allows us to estimate upper and lower error 
bounds on Hm separately, which we also report in Table 

The probability distributions also allow us to calculate 
the probability that Ha > 0 and the probability that the 
Universe is currently accelerating in its expansion for each 
SN sub-sample. The latter can be expressed via the current, 
2 = 0, deceleration parameter, go, where go = — Ha. 

As negative values of go indicate deceleration at the current 
epoch, the condition for current acceleration in the expan¬ 
sion of the Universe is given by 
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Figure 4. UPPER PANEL: The stretch-corrected SNe la Hubble diagram plotted according to the class of the host galaxy. The inset 
shows the high-redshift SNe, the main panel the Hubble diagram for the entire sample. Boxed points show high-redshift SNe excluded 
from the P99 solution. LOWER PANEL: The residuals from the adopted cosmology (‘fit-C’ of P99) for both high and low-redshift 
SNe. Excluded SNe are not plotted. Dotted-circled points discriminate between CfA SNe and C-T SNe. In both panels unfilled circles 
denote galaxies with a less secure type-classification. 
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Figure 5. Confidence regions in (Hm, S7a) for the fitting procedures in the different sub-samples listed in Table|^ The ellipses correspond 
to 68 and 90% confidence regions, and the best-fitting cosmologies are denoted by a filled diamond (general cosmology) and a filled circle 
(fiat cosmology). The diagonal line corresponds to a flat Universe, the dotted line shows the boundary between ‘no big-bang’ models, 
the dashed line denotes the infinite expansion boundary, and the dot-dashed line denotes those Universes which are currently undergoing 
acceleration. For those fits where a and Ais were fit (panels A and B), the confidence regions are obtained by integrating over these 
parameters (see text for details). 


> 




(5) 


The probabilities for these two classes of cosmological model 
are also listed in Table 

We also investigated potential environmental dependen¬ 
cies in the SN properties by plotting the Hubble diagram 
with each SN labelled according to its projected galactocen- 
tric distance (Fig. SNe located at less then 7.5kpc (the 
approximate mea n of the sample, as well as the value used 
in the analysis of [Wang et alj 1997) from the galaxy cen¬ 
tre are shown in red, those at larger distances are in blue. 
We also fit these two additional subsets of SNe using the 
same procedure as above; the results of these fits are again 
in Table ^ 

We checked the reliability of our quoted uncertainties 
in Ha (derived from the co-variance matrices of the vari¬ 
ous fits) by performing detailed Monte-Carlo simulations of 
each SN data sub-set. We created 10000 artificial datasets 
for each of the 9 sub-samples (see Table by drawing N 
SNe at a time (with replacement) from the relevant sample, 


where N is the number of SNe in each set. We then fit these 
N SNe in the same way as for the original sample. The mean 
of the resulting distribution in each fitted parameter is sat¬ 
isfactorily close to the final fitted parameter of the original 
dataset, with the standard deviation of the distribution rep¬ 
resenting approximate 1 -a uncertainties in the various fits. 
These Monte-Carlo estimated uncertainties are also listed in 
Table |^. 

We now turn to examining the significance of the vari¬ 
ous correlations in the above plots. Fig. ^ and Tablereveal 
some interesting trends. 

Firstly, each of the various samples and sub-samples 
strongly supports the SCP/HZT conclusions for the pres¬ 
ence of a cosmological constant; indeed all sub-samples ex¬ 
clude an EdS Universe at greater than the 99 per cent con¬ 
fidence level and favour cosmological models with f^A > 0 
with > 97 per cent confidence. If we restrict ourselves to 
flat flk = 0 models, this confidence increases to > 99 per 
cent for all SN sub-samples. Most noticeably, the scatter 
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Table 6. The best-fitting cosmologies for SNe residing in each of the host galaxy types. 


SNe Subset 

Median 

Best 

Best 

Monte-Carlo 

Mean 

DOF'^ 


X^ 

Probability that: 


redshift“ 



‘error’ 

dispersion^ 



/DOF'* 

Da > 0 

Da ^ 

Fits with all parameters free: 

n '79 + 0.09 
'^-0.08 








P99 ‘fit-C’ 

0.476 

0.73, 1.32 

±0.07 

0.199 (0.199) 

51 

57 

1.121 

99.79 

99.59 

P99 -1- CfA 

0.476 

0.37, 0.67 

»J-O0_o.o9 

±0.06 

0.216 (0.214) 

60 

82 

1.376 

98.24 

97.30 

Typed SNe 

0.480 

0.46, 0.77 

0.64(0.64)+°;“ 

±0.07 

0.223 (0.221) 

59 

81 

1.388 

97.96 

96.84 

Fits with a and A4b fixed: 









E/so 

0.478 

0.83, 1.02 

0.57(0.58)/:°“ 

±0.09 

0.159 (0.173) 

12 

8 

0.680 

97.09 

94.80 

Spirals 

0.480 

0.44, 0.81 

0.67{0.66)+H6 

±0.09 

0.235 (0.234) 

46 

73 

1.590 

99.85 

99.75 

Early 

0.472 

0.62, 0.65 

0.50(0.50)/;°;“ 

±0.14 

0.200 (0.193) 

22 

28 

1.309 

92.46 

86.83 

Late 

0.487 

0.45, 1.03 

0.77{0.78)/;°-p® 

±0.09 

0.272 (0.274) 

23 

40 

1.761 

99.98 

99.97 

Large Distance 

0.458 

0.72, 0.80 

f) co+0.16 
U.OC5_o 14 

±0.14 

0.196 (0.204) 

15 

19 

1.276 

94.38 

89.46 

Small Distance 

0.498 

0.77, 1.13 

^•t)0 o,o8 

±0.11 

0.256 (0.251) 

32 

56 

1.777 

99.57 

99.27 


Notes: 

“For high-redshift SNe 

^Assuming = Oi bracketed number is fit from securely classified sample 

“From best-fitting fiat cosmology; bracketed value is the dispersion from P99 ‘fit-C’ cosmology 
^For the = 0 fit* 
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Figure 6. The dispersions (in sigmas) from the best-fitting cos¬ 
mology of P99 {‘fit-C’) as a function of host galaxy type, for both 
peak magnitudes corrected for the light-curve-width-luminosity 
relation (bottom) and uncorrected magnitudes (top). The colour 
coding is as in Fig. 


consider the securely classified sample of host galaxies (see 
Table P) . We will return to the significance of this conclusion 
in the next section. 

The reduced scatter we observe for the E/SO sample is 
an important result. Accordingly, we are interested in under¬ 
standing whether it is a robust result. To verify this, we have 
performed simulations to determine how likely we would be 
to obtain similar sub-samples to the E/SO sample by chance. 
We draw Veil high-redshift SNe from the full sample (where 
Veil is the number of high-redshift SNe in E/SO galaxies) 
and perform the same best-fitting procedure to these SNe. 
By repeating the test a large number (10000) of times, we 
investigated how frequently the derived cosmological param¬ 
eters produce a Ha smaller than the actual E/SO sample and 
a smaller mean dispersion about the best-fitting cosmology. 

We find that less than 3 per cent of the random samples 
possess a Qa smaller than that found for E/SO galaxies, and 
a mean dispersion about the best-fitting cosmology smaller 
than that for the E/SO SNe. Thus it seems the E/SO sample 
is unlikely to be simply a sub-sample selected randomly from 
the larger sample of SNe, and that there is a significant 
difference in the photometric properties of the E/SO SNe 
that differs from the mean of the full SNe sample. 

Finally, we can see no clear trends in the Hubble dia¬ 
gram plotted according to the projected distance of the SN 
from the host galaxy (Fig. ^), again regardless of whether 
stretch corrections are applied. 


about the best-fitting line for SNe arising in spiral or irreg¬ 
ular hosts is greater than that for SNe arising in spheroidal 
host galaxies (see also the x^/DOF in Table This trend 
of increased scatter with class is particularly striking and 
can be seen more clearly by examining the dispersions from 
the best-fitting cosmology (Fig. 

If patchy dust extinction were the origin of the increased 
scatter in the later type hosts, we would expect, on aver¬ 
age, the SNe in these systems to appear fainter. Such a 
trend is indeed seen; SNe occurring in spiral galaxies ap¬ 
pear fainter, with or without the application of a stretch- 
correction. These conclusions are not affected if we only 


5.2 Correlations between light-curve parameters 
and environment 

In st udies of the environments of local SN e la, many authors 
(e.g. Riess et al. 1999; Hamuy et al. 200C) find that intrinsi¬ 


cally brighter SNe (SNe with larger stretch values) occurred 
in late-type galaxies, and intrinsically fainter (small stretch) 
SNe were preferentially located in spheroidal hosts. It is im¬ 
portant to ensure that the effects we see in Fig. ^ are not 
systematic effects simply related to the application stretch 
corrections which vary accordi ng to the type of the host 
galaxy (e.g. Hamuy et al. 2000). 
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Figure 7. As Fig. but with the SNe plotted according to the projected location from the centre of the host galaxy. Stretch corrections 
are applied. Filled circles show SNe lying at < 7.5 kpc, whilst empty circles show SNe at larger projected distances. 


We have considered this issue in two wa ys. First, we 
have repeated the cosmological fits of Section ^d|, compar¬ 
ing the dispersions from the best-fitting cosmology with and 
without applying the stretch corrections to each peak lumi¬ 
nosity (Fig. We find that the trends in Fig. ^ remain 
even if stretch corrections are not applied, with very simi¬ 
lar derived cosmological parameters within the uncertainties 
quoted in Table 


More directly, we can examine the distribution of the 
stretch parameter for each SN located in each galaxy type. 
Fig. 1^ shows this compariso n for the hig h and low redshift 
sam ples. Th e trend noted by Riess et al (1999) and Hamuy 
et al. ( 1000) is apparent in the low-redshift sample: an ab¬ 
sence of small-stretch SNe in late-type galaxies and of large- 
stretch SNe in early-type galaxies . However, this same trend 
is not seen in the SCP high-redshift SNe sample. This is due 
to the fact that, as shown in P99 (their fig. 4), the stretch 
distribution of the high-redshift SNe is narrower than that 
seen in the low redshift samples. We return to this in Sec¬ 
tion ^ 

Finally, Fig. ^ shows the SN stretch plotted as a func¬ 
tion of projected distance from the host galaxy. Here a weak 
trend is found possibly suggesting that SNe at greater pro¬ 
jected distance s may be fa i nter, similar to that seen in the 
local sample of jRiess et"^ ( 1999 ), This could possibly arise 
if, within a distribution of luminosities, fainter SNe were 



6 DISCUSSION 

Of the various results presented in Section the most robust 
and significant trend we have found is that SNe found in spi¬ 
ral and irregular galaxies present a larger scatter around the 
best-fitting cosmological model than those located in E/SO 
galaxies. However, SNe in all galaxy types provide evidence 
for a significant cosmological constant (at ~ 5a). Addition¬ 
ally, SNe in late-type galaxies also appear fainter than those 
found in E/SOs, regardless of whether stretch-corrections are 
applied. 

An obvious explanation of these main results, which do 
not depend on whether stretch corrections are made, is that 
there is increased extinction in later type galaxies and of 
an amount which varies somewhat from galaxy to galaxy. 
Though the SNe magnitudes are corrected for Galactic ex¬ 
tinction, no account is taken of dust residing either in the 
host galaxy or along the line of sight. 

P99 investigated the effect of differential reddening on 
the high-redshift SNe as compared to their local counter¬ 
part by contrasting the mean B — V colour of high redshift 
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Figure 8. The distribution of SN stret ch with host galaxy type. 
The low-red s hift S Ne are taken from the Hamuy et al. (1996a) and 
(|l999 ) samples; see text for further details. All SNe 


Riess et al. 


(including those rejected as stretch outliers from the cosmological 
fits) are shown in this figure. 



Projected Distance (kpc) 


Figure 9. The SN stretch plotted as a function of projected 
distance from the centre of the host galaxy. The colour coding is 
as in Fig. ^ 


SNe at maximum light with that for the local sample. They 
find that the colour excess distributions show no significant 
differences at low and high redshift, with the reddening dis¬ 
tributions for the high-redshift objects consistent with the 
reddening distribution for low-redshift SNe within the mea¬ 
surement uncertainties, and both samples possessing almost 
identical error-weighted means. Additionally, they perform 
cosmological fits with each SN corrected for extinction ac¬ 


cording to its E(B — V), and find that the fitted parameters 
change little with this modification. 

In the following we have revisited and extended this 
analysis. Although the precision is poor for an individual 
SN, it is advantageous to re-examine the situation as we now 
have independent diagnostics which relate to the presence 
of host galaxy dust. The differential reddening affecting any 
one SN can be determined by comparing its observed Bmax — 
Knax with the (Bmax — 14iax)o for low-redshift SNe with the 
same stretch-factor. The differential colour excess E{B — V) 
is determined as 

E{B -V) = (Bmax - Knax)obs “ (Bmax “ Knax)o- (6) 


The (fJmax — V'max)obs are calculated following P99. For 
SNe in the redshift range 0.3 < z < 0.7, the R and I band 
measures are used instead of B and V and (B max — Vmax)obs 


is calculated using a 

ppropriate /i:-corrections flKim, Goobar, 

& Perlmuttei 

19961 

Nugent, Kim, & Perlmutter 2002). For 


SNe at j* > 0.7 the R and I measures correspond more 
closely to SN rest-frame U and B bands, so E{B — V) is 
calculated from E{U — B); likewise for SNe at z < 0.3, where 
the R and I correspond to rest-frame V and R, E{B — V) is 
calculated from E{V — R). We omit those 6 of the 42 high- 
redshift SNe that do not have both R and I measures from 
this analysis. 

The distribution of this E{B — V) parameter as a func¬ 
tion of the host galaxy type is shown in Fig. p^ . The errors 
in E{B — V) are derived from the uncertainties in the R and 
1 photometry. The error-weighted mean values of E{B — V) 
for each class are -0.01, 0.20 and 0.03 for E/SO, early spirals 
and later spirals/irregulars respectively; the value for spirals 
taken as a whole is 0.10. We can see that SNe lying in E/SO 
galaxies form a very narrow and tight distribution around a 
mean E{B — E) ~ 0, whereas SNe residing in spiral galax¬ 
ies show a larger range of E{B — V) values and means of 
E{B — E) > 0. This is suggestive that, as expected, SNe in 
E/SO galaxies are unlikely to be affected by the presence of 
dust, whereas the effect of dust on SNe in spiral systems is 
more uncertain. Eurthermore, the two SNe rejected from ‘fit- 
C’ in P99 as ‘likely-reddened’ (1996cg and 1996cn) are both 
located in early spiral galaxies; indeed, these two objects are 
the primary cause of the large mean value of E{B — V) found 
for these systems. 

The above result is, to some extent, independent of 
the morphological trends seen in the Hubble diagram since 
the former relies on the SN colour whereas the latter de¬ 
pends on the peak luminosity. Naturally there are possi¬ 
ble inter-dependencies but we argue both results provide 
strong support for our two main conclusions. Eirstly, much 
of the scatter in the SNe la Hubble diagram arises from ex¬ 
tinction in the later type galaxies. As this scatter is host 
galaxy morphology-dependent, clearly it is internal to the 
host galaxies. Secondly, not only do we see a tight scatter 
in the Hubble diagram for SNe in E/SOs, but differentially 
there is no obvious difference in the extinction as inferred 
from the SNe at high and low redshift. Accordingly, the cos¬ 
mological conclusions derived from this E/SO subset must be 
particularly robust to systematic errors arising from extinc¬ 
tion. As can be seen in Table this data provides a ~ Sa 
detection of a non-zero cosmological constant. 

An interesting question then is; how much extinction 
in high-redshift SNe la is implied by our results? We can 
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Table 7. The best-fitting AA b for SN la arising from different galaxy types, assuming the best-fitting cosmology for the E/SO sub-sample. 



E/so 

Early Spi 

Late Spi/Irr 

All Spirals 

Mb 

-3.39 ±0.07 

-3.37 ±0.08 

-3.25 ±0.06 

-3.32 ±0.05 

Difference with E/SO sample 

- 

0.02 ± 0.11 

0.14 ±0.09 

0.07 ±0.09 



Figure 10. The SN rest-frame ‘colour excess’ {E{B — V)) plotted 
as a function of the host galaxy type. E{B — V) is as defined in the 
text (see also P99). The error-weighted mean of each distribution 
are shown as stars; the error bars are derived from the R and I 
photometry for each SN. The vertical offsets for each plotted SN 
is used to avoid over-plotting the data points and is illustrative 
only. Only the high-redshift SNe are shown. 


address this question in the following manner. We repeat 
our cosmological fits for the E/SO and various spiral SNe 
sub-samples, but this time fit for the parameter AAb rather 
than flM and Qa, and instead hold Hm and SIa fixed at their 
values as derived from the (presumed dust-free) E/SO sample 
of SNe. As Ads is the intrinsic B-band peak magnitude of a 
SN la, if those SNe in E/SO galaxies are less affected by dust 
than those in spiral galaxies, we might expect the best-fitting 
Ads to be brighter in these galaxies. These best-fitting Ads 
values and l-cr errors are given in Table 

Assuming the best-fitting E/SO flat-Universe cosmol¬ 
ogy, SNe appear 0.14 ± 0.09 mag fainter in late-type spirals 
than in E/SOs. However, these uncertainties are currently 
large (due to the small number of objects in the different 
samples), and the result is consistent with both a near-zero 
offset and an offset of 0.23 mag. Though there are obvious 
caveats associated with this result, and the scatter from SN 
to SN is large due to the photometric uncertainties in the SN 
magnitudes, this measure would however appear to suggest 
that the amount of extinction suffered by these high-redshift 
objects is small, confirming earlier conclusions (e.g. P99). 

This result a ppear s in contradiction to the work of 
Rowan-Robinson (p002|) , who suggest that the evidence for 
a positive cosmological constant is reduced if an average 
internal dust-extinction correction of 0.33 mag for all high- 


redshift SNe is made. However, median extinction values as 
large as this do not appear consistent with these new results, 
particularly in the E/SO systems where internal extinction 
is expected to be very small, and which show that Ha > 0 
at ~ 5cr. 

Modest extinction for even the SNe occurring in late- 
type galaxies is supported by the small amount of spectral 
data available for the emission line host galaxies (Table 
though the extinctions implied are not as small as those in¬ 
dicated by the dispersions in the Hubble diagram. Eor the 6 
galaxies where a measure of the Balmer decrement is possi¬ 
ble (all at z < 0.5 due to Ha shifted out of the ESI spectral 
coverage) the measures suggest (except for one case) small 
extinctions on the nebular emission lines, and hence smaller 
extinction on the stellar continuum l ikely affecting t he ob - 
served SNe, using the prescription of Calzetti et al. (2000). 
However, most of the measured Av are still larger than the 
0.14 mag implied by the difference in Mb between E/SO and 
spiral systems (Table We hypothesise this is due to SNe 
lying well away from the H ll regions from which the nebular 
lines used to measure the dust content originate, an effect 


which is not removed even after using the Calzetti et al 
( poo^ prescription. In the one case where the extinction 
measured is large {Ay — 3, 1995az), the SN was located at 
a large projected galactocentric distance of ~ 13 kpc, one 
of the largest in the sample (see Fig.s and where the 
large amounts of dust that were detected in the H ll regions 
may not be present. Clearly larger samples of spectra are 
required to confirm these trends. 

At first sight, the small extinction found for SNe in late- 
type galaxies appears a surprising result. Observational evi¬ 
dence suggests that the amount of star-fo rmation in norma l 
field galaxies increase with redshift (e.g Lilly et al. 1996), 


and that much of this increase results fro m late-type and 
irregular systems (Brinchmann et al. 1998), which presum¬ 


ably should be dustier systems than their non-starforming 
counterparts, as increased star-formation activity is usually 
associated with increased dust content. Although we see 
such a correlation in our data, in absolute terms the SNe 
of P99 even in the later-type hosts appear to suffer little 
extinction, as confirmed by the limited spectroscopic data 
presented here and the colours of the low and high-redshift 
SNe at maximum light. 


Hatano et al. (1998) have investigated the effects of ex¬ 
tinction on the optical properties of SNe. For simple assump¬ 
tions concerning the spatial distribution of dust and SN pro¬ 
genitors in a model disc galaxy, they find that the amount 
of extinction suffered by SNe is small. The typical ranges in 
mean Ab is 0.3 mag to > 1 mag when all SNe are considered; 
however these means are dominated by a few high extinction 
events - most SNe are mildly obscured, but with a long tail 
to higher extinction events, which drop out of flux-limited 
surveys such as the SCP/HZT search campaigns. Mean val¬ 
ues for an ‘extinction-limited subset’ (Ab < 0.60) range 
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from 0.12 to 0.16 mag, depending on the inclination of the 
system, in good agreement with the 0.14 mag difference we 
hnd between E/SO and late-type spiral SNe. 

A further explanation for the low extinction on high- 
redshift SNe in late-type systems is that the progenitors of 
SNe la are largely located in the outermost components of 
disc and irregular galaxies, perhaps as landmarks of earlier 
stages of galactic assembly, whereas the dust is distributed 
within regions of more recent star-formation activity con- 
hned to the thin disc. In popular hierarchical models, where 
the bulk of galactic assembly occurs over 1 < z < 3, one 
would expect much mixing of dust and old stars. The ab¬ 
sence of signihcant extinction in the bulk of our sample 
might be used to argue against a significant increase in the 
merger rate prior to a ~ 0.5 — 0.8. 

Another possible explanation is related to the selection 
effects inherent in SNIa surveys; do SNIa searches sample 
the entire range of galaxy types available at a given red- 
shift? Alternatively, are we missing SNe that reside in par¬ 
ticularly dusty host galaxies? In principle, we can examine 
this by comparing the morphological mix of our host galaxy 
sample to the observed mix of normal galaxies in the same 
redshift range. The most relevant study to date is that of 


Brinchmann et al. 


(1998), who analyse HST imaging of the 


galaxies selected in the CFRS and LDSS redshift surveys. 
They classify their sample into three broad classes - E/SOs, 
spirals and irregulars - and find a large increase with red¬ 
shift in the fractional number of irregular objects from 9 per 
cent at z < 0.5 to over 30 per cent at 0.7 < 2 < 0.9, with de¬ 
creases in the relative numbers of E/SOs and regular spirals. 
This leads to a change in t he split of both integra ted B- 
light luminosity density (see Brinchmann et al. 1998 ), and a 
change in th e morphological split of in tegrated stellar mass 
density (see Brinchmann & Ellit 200C), both likely related 
to SNIa events. Hence, in a truly volume-limited survey of 
SNe la, we might expect the fraction of events found in ir¬ 
regular galaxies to sharply increase at 2 > 0.5, tracing the 
increased stellar mass in these systems. 

Comparisons with this study are obviously problematic 
due to the small size of the current sample and the differing 
redshift ranges, as well as differing classihcation techniques 
for the galaxies. We have 16 securely classihed galaxies at 
0.5 < 2 < 0.8, of which 7 are classed as type-2 (44 ± 17 per 
cent). However, the Brinchmann et al. scheme would class 
many of these regular late-type galaxies separately from 
those which are irregular systems - only 3 of the objects 
(19 ± 11 per cent) are irregular in this definition scheme. 
Brinchmann et al. find a fraction of 33 ±6 per cent, opening 
the possibility that we do not see this same rapid rise in the 
number of ‘truly’ irregular systems in the SN sample. 

This possible lack of irregular systems can be recon¬ 
ciled if we hypothesise that such objects contain, on average, 
larger amounts of dust due to their presumably higher star- 
formation activity. This would result in any SNe located in 
these galaxies to be dimmed to a greater extent than those 
in E/SO or regular spiral systems; the SNe searches could 
therefore be (implicitly) biased against extincted SNe. Con¬ 
ceivably, many SNe la do suffer signihcant extinction but, 
because of the flux-limited SNe search criteria, they are not 
included in our sample. 

Our observation that SNe la residing in E/SO galaxies 
show a smaller scatter about the best-fltting cosmological 


model (as revealed by the x^/DOF for the various flts in 
Table ^ allows us to estimate the intrinsic spread in the 
peak magnitude of these SNe compared to those residing in 
other classes of galaxy. Currently, each SNe error contains 
a contribution of 0.17 mag added in quadrature to allow for 
the uncer tainty in the light-c urve-width-luminosity correc¬ 
tion (e.g. Hamuy et al 1996b), regardless of the host class. 


As some of this dispersion arises from host galaxy extinction, 
we can use the (presumed) dust-free E/SO SNe population 
to estimate the intrinsic dispersion for this class of SNe. For 
SNe in each host type, we therefore repeat the cosmological 
fits, altering the intrinsic dispersion until the x^/DOF for 
each fit is equal to one. We perform this test for the high 
and low-redshift SNe taken together, and also solely for the 
high-redshift SNe. 

The intrinsic dispersions are (in magnitudes) 0.09(0.10), 
0.29(0.25), and 0.69(0.24) for E/SO, early and late-type spi¬ 
rals/irregulars respectively - bracketed values refer to fits 
where only high-redshift SNe are used (the values for spirals 
taken as a whole are 0.46 and 0.26) . Note that the small 
intrinsic dispersion seen for those SNe in E/SO matches very 
closely the dispersion seen in low-redshift SNe after a tech¬ 
nique correcti ng for the effects of r eddening has been applied 
(0.11 mag; see Phillips et al 1999), and the large dispersion 


seen in late-type spiral classes is heavily influenced by one 
or two low-redshift outliers. The implication from this small 
dataset is clear: SNe residing in E/SO hosts are superior 
‘standard candles’ than those in later-type galaxies due to 
the small amount of extinction affecting them. 

Finally, it is important not to overlook a possibly im¬ 
portant result suggested by the lack of range in light-curve 
width (or stretch) indicated in Fig. ^ This hgure does pro¬ 
vide some evidence of a shift in the population distribution 
between SNe found at low and high-redshift: at high-redshift 
we see no trend in stretch with galaxy type, whereas at low- 
redshift SNe with larger stretches are found in later-type 
galaxies and are missing from E/SO galaxies. It’s impor¬ 
tant to emphasise that this population shift has no impact 
on the security of the cosmological results, as light-curve- 
width differences are accounted for in the stretch correction 
of equation (|^, and stretch outliers are in any case rejected 
from the flts. Given the small number of SNe in each group. 
Fig. ^ will, of course, be sensitive to one or two extreme or 
unusual events. Bearing this caveat in mind, we performed 
Kolmogorov-Smirnov (K-S) tests on the different stretch 
distributions for SNe residing in each galaxy class, and And 
the significance levels for the hypothesis that the SNe at low 
and high-redshifts are drawn from different distributions to 
be 89, 71 and 7 per cent for E/SO, early- and late-type galax¬ 
ies respectively (the significance level for low-redshift SNe in 
E/SO and late-type galaxies having a different stretch distri¬ 
bution is 94 per cent in this sample). Though these trends 
are only suggestive given the small number of events, they 
might be explained if we assume that the light-curve proper¬ 
ties of SNe la depend primarily on the age of the progenitor 
system. 

At low redshift, it appears intrinsically fainter SNe are 
located in older (E/SO) stellar environmen ts rather than 
young spiral systems (e.g. Hamuy et al. 2000 ). As we look to 
higher-redshift (i.e. viewing galaxies at an earlier stage in the 
evolution of their composite stellar populations), the E/SO 
galaxies will posses younger populations than is the case in 
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the local Universe (at 2 = 0.5, the age of the Universe is 
~ 8 Gyr compared to ~ 13 Gyr at 2 = 0 in a AGDM cos¬ 
mology) . If, as was discussed in Section ^ older progenitors 
generate fainter SNe events, at high-redshift we would not 
expect to see the fain ter (smaller stretch ) SNe found in local, 
older, E/SO galaxies (Hamuy et al. 200C), thus strengthening 
the case for a cosmological constant. 


7 CONCLUSIONS 


We have investigated the Hubble diagram of Type la super¬ 
novae spanning a large range in redshift and have classified 
39 distant events according to the nature of their host galax¬ 
ies as revealed by HST imaging, intermediate dispersion 
spectroscopy and broad band colour. Together with mor¬ 
phological data for the host galaxies of 25 local SNe drawn 
from the literature, we find the following. 

(i) The scatter on the SNIa Hubble diagram, as deter¬ 
mined from residuals of the peak brightness with respect to 
the best fit cosmological model, correlates closely with the 
host galaxy type. The scatter is minimal for SNe occurring 
in host galaxies broadly classed as early-type, and increases 
towards later classes. The correlation does not depend on 
whether the SN are corrected for differences in light-curve 
shape and thus cannot arise primarily from type-dependent 
light-curves. 

(ii) The above trend is accompanied by the observa¬ 
tion that SNe occurring in the later types are on average 
marginally fainter than those in E/SO galaxies, with a dif¬ 
ference in their absolute peak lumi nosities of 0.14T0. 0 9 mag . 

(iii) Following the approach of Perlmutter et al, (1999), 
we compare the colour-excesses of distant SNe la at maxi¬ 
mum light with those observed locally for examples of the 
same stretch and see an indication of a type-dependent 
colour excess of a similar nature, with, on average, spirals 
hosting redder SNe. 

The data are also suggestive of a lack of any significant 
correlations at high-redshift between the light-curve stretch 
parameter, which correlates with the observed peak luminos¬ 
ity, and the host galaxy type. The full implications of this 
result are currently unclear, given the (current) small sample 
sizes. However, the lack of small stretch (fainter) SNe and 
the increase in larger-stretch SNe in high-redshift E/SOs, 
compared to that found in low-redshift samples, suggests 
that the age of the progenitor system may be the primary 
variable in determining the light-curve properties of a par¬ 
ticular SN, with younger systems located at higher-redshift 
when compared to those seen locally. It is important to note 
that these potential trends will not affect the validity of the 
cosmological result. Larger samples of SNe are clearly re¬ 
quired to probe these trends further. 

We find weak trends between light-curve shape and the 
spatial location of the SNe in the host galaxy, as claimed oc¬ 
casionally in lower redshift samples, in the sense that fainter 
SNe are located away from galaxy centres. However, there 
are two caveats associated with this result: firstly, there is 
a selection effect of faint SNe being missed near to galaxy 
cores at low-redshift (this effect is not beli eved to operate t o 
the same extent at high-redshift, see e.g. Pain et al. 1996), 
and secondly, we use projected distances rather than de- 
projected distances in the analysis. 


We interpret the above results according to the hypoth¬ 
esis that dust extinction in the host galaxy is the significant 
cause of the scatter in the Hubble diagram. The implications 
of this conclusion are significant. Most importantly, the Hub¬ 
ble diagram eonfined to SNe la oceurring in early-type host 
galaxies with presumed minimal internal extinction presents 
a very tight relationship and provides a ~ 5 ct confirmation of 
a non-zero cosmological constant, assuming a flat Universe. 
When no assumption about the flatness of the Universe are 
made, we find that SNe in E/SO galaxies imply Ha > 0 at 
nearly 98 per cent probability. 

The differential trends in extinction by host galaxy mor¬ 
phology as revealed by the SN colour excess and mean offset 
in the type-dependent Hubble diagram with respect to that 
defined by the early-type hosts implies only modest extinc¬ 
tions in late-type spirals of Ab = 0.14 ± 0.09 mag, though 
with a large galaxy-to-galaxy scatter. We discuss briefly the 
origin of this small extinction in the light of the increasing 
fraction of star-forming and presumed dusty galaxies at high 
redshift and suggest a number of possible explanations. 

Our study is a good illustration of how ancillary data on 
host galaxies can be used to examine the nature and reliabil¬ 
ity of SNe as probes of the cosmological expansion history. 
Gontrary to suggestion s made by some previous studies (e.g. 


Rowan-Robinson 


2002), we show that internal dust extinc¬ 


tion cannot be a primary contaminant and that supernovae 
of Type la, particularly those occurring in spheroidal galax¬ 
ies which can be readily screened with spectroscopic, colour 
and morphological data, represent a very powerful cosmolog¬ 
ical probe. This conclusion is of particular significance for 
future projects such as the SuperNova Acceleration Probe 
(SNAP) experiment (e.g. Aldering et al. 2002), a space- 
based telescope designed to determine the cosmic expansion 
history to high precision via studies of many thousands of 
SNe la. This instrument will provide HST-quality imaging 
and nine-band optical and near-infrared colour photometry 
for the host galaxies of nearly every SN it discovers, allowing 
this host-galaxy study to be extended using better statistics 
and a superior classification of the host-galaxy types, and a 
particularly robust cosmological solution from those events 
located in spheroidal systems. 
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